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Constraining the geometry of the neutron star RX J1856.5-3754 
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ABSTRACT 

RX J1856. 5-3754 is one of the brightest, nearby isolated neutron stars, and considerable ob- 
servational resources have been devoted to its study. In previous work, we found that our latest 
models of a magnetic, hydrogen atmosphere matches well the entire spectrum, from X-rays 
to optical (with best-fitting neutron star radius R x 14 km, gravitational redshift Zg ~ 0.2, 
and magnetic field Bs!4xlO'~G). A remaining puzzle is the non-detection of rotational 
modulation of the X-ray emission, despite extensive searches. The situation changed recently 
with XMM-Newton observations that uncovered 7 s pulsations at the ^ 1 % level. By compar- 
ing the predictions of our model (which includes simple dipolar-like surface distributions of 
magnetic field and temperature) with the observed brightness variations, we are able to con- 
strain the geometry of RX J1856. 5-3754, with one angle < 6° and the other angle a; 20-45°, 
though the solutions are not definitive given the observational and model uncertainties. These 
angles indicate a close alignment between the rotation and magnetic axes or between the ro- 
tation axis and the observer. We discuss our results in the context of RX J 1856.5-3754 being 
a normal radio pulsar and a candidate for observation by future X-ray polarization missions 
such as Constellation-X ox XEUS. 

Key words: polarization - stars: individual (RX J1856. 5-3754) - stars: magnetic fields - 
stars: neutron - stars: rotation - X-rays: stars 



1 INTRODUCTION 

Seven candidate isolated, cooling neutron stars (INSs) have been 
identified by the ROSAT All-Sky Survey (see Kaspi, Roberts, & 
Harding 2006; Haberl 2007; van Kerkwijk & Kaplan 2007, for 
recent reviews). These objects share the following properties; (1) 
high X-ray to optical flux ratios of log(/x//opticai) ~ 4 - 5.5, (2) 
soft X-ray spectra that are well described by blackbodies with 
kT ~ 50 - 100 eV, (3) relatively steady X-ray flux over long 
timescales, and (4) lack of radio pulsations. 

For the particular INS RX 11856.5-3754, X-ray and opti- 
cal/UV data can be well-fit by two corresponding blackbody spec- 
tra; the X-ray spectrum by a blackbody with kT^ = 63 eV 
and emission size = 5.1(d/140pc) km, and the optical/UV 
spectrum by a blackbody with kT^f = 26 eV, and R^^f = 
21.2 (rf/ 140 pc) km (Burwitz et al. 2001; van Kerkwijk & Kulka- 
rni 2001a; Drake et al. 2002; see also Pavlov, Zavlin, & San- 
wal 2002; Trumper et al. 2004; Walter 2004), where = refr/(l + 
Zg), R'^ = R™(1 + Zg), R™ is the physical size of the emission re- 
gion, and d is the distance. The gravitational redshift Zg is given by 
(1 -I- Zg) = (1 - 2GM/Rc^)-^'^, where M and R are the mass and 
radius of the NS, respectively. The high temperature, small area X- 
ray blackbody suggests a small hot spot on the NS surface, while 
the optical blackbody with R^^^ is the remaining large, cool surface. 

Even though blackbody spectra fit the data, one expects NSs 
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to possess atmospheres of either heavy elements (due to debris 
from the progenitor) or light elements (due to gravitational set- 
tling or accretion). The lack of any significant spectral features 
in the X-ray spectrum argues against a heavy element atmosphere 
(Burwitz et al. 2001, 2003), whereas non-magnetic or fully-ionized 
magnetic hydrogen atmosphere spectra do not provide a good fit 
(Pavlov et al. 1996; Pons et al. 2002; Burwitz et al. 2001, 2003). 
For the last case, it is important to note that, since kT ~ tens of eV 
for RX J1856.5-3754 and the ionization energy of hydrogen at 
B = 10'^ G is 160 eV, the presence of neutral atoms must be 
accounted for in the magnetic hydrogen atmosphere models; the 
opacities are sufficiently difi'erent from the fully ionized opacities 
that they can change the atmosphere structure and continuum flux 
(Ho et al. 2003; Potekhin et al. 2004), which can afi'ect fitting of the 
observed spectra (see, e.g.. Ho et al. 2007). 

Another issue that may argue against the two-temperature 
blackbody model is the the non-detection, until recently, of X-ray 
pulsations due to the rotation of the NS (down to the 1.3% level; 
Drake et al. 2002; Ransom, Gaensler, & Slane 2002; Burwitz et 
al. 2003; Zavlin 2007). As discussed in the important work by Braje 
& Romani (2002), very low rotational pulsations could be due to a 
very close alignment between the rotation and magnetic axes and/or 
an unfavorable viewing geometry. Assuming blackbody emission, 
they find this is possible in ~ 5% of viewing geometries (with 
R = 14 km). However, more realistic atmosphere emission can lead 
to enhancement of pulsations, as well as an energy-dependence of 
the pulsations (see below). 
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Thus there existed three major observational and theoretical 
inconsistencies: (1) the inferred emission size from blackbody fits 
are either much smaller {R°° ~ 5 km from the X-ray data) or much 
larger ~ 20 km from the optical/UV data) than the canonical 
NS radius of 10 - 12 km, even after correcting for the redshift, (2) 
blackbodies fit the spectrum much better than realistic atmosphere 
models, and (3) strong upper limits on X-ray pulsations suggest 
RX J1856.5-3754 may have a largely uniform temperature over 
the entire NS surface. In response to these problems, we applied 
our latest magnetic, partially-ionized hydrogen atmosphere mod- 
els and obtained a good fit to the entire multi-wavelength spectrum 
of RX J1856.5-3754 (Ho et al. 2007; hereafter H07). Besides us- 
ing the more realistic atmosphere, the fit requires a smaller emis- 
sion size R^™ X 14 km (as compared to the blackbody models); 
if interpreted as emission from the entire NS, the size can be sat- 
isfied by a stiff but standard equation of state (see, e.g., Lattimer 
& Prakash 2007). Furthermore, we found geometries that produce 
pulsations near the observed limits. Because of observational un- 
certainties and the computationally tedious task of constructing a 
complete grid of models, the results may not be unique. Never- 
theless, our model represents the most self-consistent picture for 
explaining all the observations of RX J1856.5-3754. 

Very recently, XMM-Newton observations uncovered pulsa- 
tions from RX J1856.5-3754 with a period of 7 s (pulse amplitude 
K 1%) and an upper limit on the period derivative; by assuming 
vacuum magnetic dipole braking, this implies B < 10'"* G (Tiengo 
& Mereghetti 2007; hereafter TM07). The detection of pulsations 
allows us to constrain the geometry (i.e., two angles a and ^, where 
a is the angle between the rotation and magnetic axes and ( is the 
angle between the rotation axis and the direction to the observer) 
through a comparison of the observed brightness variations (= light 
curves) with those implied by the model of H07; the latter are ob- 
tained by computing the (rotation) phase-dependent spectra from 
the entire NS surface. 

Spectra from the whole NS surface are necessarily model- 
dependent (see, e.g., Zavlin et al. 1995; Zane et al. 2001; Ho & 
Lai 2004; Zane & TuroUa 2006), as the magnetic field and tempera- 
ture distributions over the NS surface are unknown. Magnetic field 
variations over the surface will induce surface temperature vari- 
ations (Greenstein & Hartke 1983). Furthermore, radiation from 
the surface of a magnetic NS with an atmosphere differs signif- 
icantly from that implied by isotropic blackbody emission: the 
emission depends on the direction of the local surface magnetic 
field (Shibanov et al. 1992; Pavlov et al. 1994). Indeed, there have 
been recent works attempting to fit magnetic atmosphere spectra 
to observations of NSs other than RX J1856.5-3754. Lloyd et 
al. (2003) fit the spectrum and light curve of PSR 0656-1-14, while 
Zavlin & Pavlov (2004), using a different model, do the same for 
PSR B0950+08. Zane & Turolla (2006), using the method of Lloyd 
et al. (2003) to produce an extensive library of light curves for fully 
ionized hydrogen atmospheres, fit the light curves of several INSs. 

An outline of the paper is as follows. In Section[2l we describe 
the model used to produce light curves of RX J1856. 5-3754. Light 
curves for arbitrary angles (a,^) are discussed in Section [3] The 
observations and fits to the observations are shown in Section |4] 
We summarize and discuss our results in Section[5] 

2 MODEL FOR NEUTRON STAR EMISSION 

We adopt a relatively simple model (see also Zavlin et al. 1995; 
Lloyd et al. 2003) for the surface magnetic field B and effective 



Table 1. Parameters for Neutron Star Surface 
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temperature T^ff distributions: we assume the surface is symmet- 
ric (in B and T^if ) about the magnetic equator and divide the hemi- 
sphere into four magnetic colatitudinal regions. We generate (local) 
atmosphere models (see H07, and references therein, for details; 
see Zavlin & Pavlov 2002, for a review of atmosphere modeling) 
for each region with the parameters given in Table [T] where 0^ 
is the angle between the local magnetic field and the surface nor- 
mal. Emission from any point within a colatitudinal region is given 
by the atmosphere model for that region. Note that the magnetic 
field distribution is roughly dipolar. The phase-resolved spectra and 
light curves from the entire NS surface is then calculated by the 
method described below (see the analogous formalism in Pavlov & 
Zavhn 2000). We assume M = IAMq and = 14 km (see H07). 
This implies Zg = 0.2, and all energies quoted in this work have 
been redshifted from the local NS frame by this factor. 

The observed emission from a local region (e.g., the hot mag- 
netic polar cap) of a rotating NS depends on two angles: the angle 
a between the rotation and magnetic pole axes and the angle ( be- 
tween the rotation axis and the line-of-sight to the observer. In the 
coordinate frame x, where z is the direction to the observer and the 
spin axis il, is in the xz-plane, the polar angle G„, and azimuthal 
angle (f),„ of the magnetic pole m are given by 

cos9„, = cos 0p sin cf sin f + cos Q- cos ^ (I) 
sin 0p sin a 

tan<p,„ = , (2) 

cos c/>p sin a cos ( - cos a sin ^ 

respectively (see Figure [T). 0p (= Q.J) is the rotation phase. Note 
that the light curves are symmetric under the interchange of a and 
since the brightness depends only on the angle d,„ between the 
magnetic and line-of-sight axes. At the emission-point P [=Pi9, 0)] 
on the NS surface (see Figure[T), the direction of the local magnetic 
field B (assuming a magnetic dipole geometry) is given by 

cosQg = cos(0,„ - 0) sin 6 sin S,,, -I- cos cost/,,, (3) 
^ , sin(0,„ -0)sin6',„ 

tan^i = . (4) 

cos(0„i - (p) cos 6 sin 6,,, - sin 6 cos 6,,, 

(f>i; is the azimuthal angle between B and the photon wavevector k 
and is defined in the coordinate frame x„, where z,, is the direction 
of the surface normal and k is in the XnZn-plane. 

To determine the flux spectrum Fe for a given rotation phase, 
we first compile a table of specific intensities lE(Sk,<f'k), where 9^ 
is the angle between k and z„. are computed from atmosphere 
models described in H07; these models depend on the local effec- 
tive temperature and magnetic field direction &b and strength (see 
Table[T). The flux spectrum can then be calculated from 

^^ = (t) (r4j^I '^^^I (5) 

where yu^ = cos dk. In the absence of the bending of the path of light 
due to gravity, the polar angles of the emission-point and the photon 
wavevector are the same, i.e., 6 = 9^. Gravitational light-bending 
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Figure 1. Coordinate axes and angles used to describe tlie pulsed emission. 
Tlie observer is in tlie z-direction. The dashed line indicates the position of 
the magnetic pole m as the neutron star rotates about £is- 



causes more of the NS surface to be visible (for M = IAMq and 
R = 14 km, 6 < 115° as compared to < 90° without light-bending). 
We use the approximate relation between 6 and 6^ given in Be- 
loborodov (2002) 



cos 6 - 



cos - IGMjc^R 



1 



(6) 



IGM/c'-R ' 

which deviates from the exact relation given in Pechenick, Ftaclas, 
& Cohen (1983) by < 1% for our chosen M and R. The light curves 
are then computed by integrating the photon count spectra over 
the energy range of the XMM-Newton observations (see TM07): 
0.15-1.2 keV for the total X-ray spectrum, 0.15-0.26 keV for the 
soft band, and 0.26-1.2 keV for the hard band. 



3 RESULTS 

Figure |2] shows the light curves of our NS model for various ge- 
ometries (cf, f). We also plot the analytic light curves from Be- 
loborodov (2002) for isotropic emission from two antipodal hot 
spots (see Zavlin & Pavlov 1998; Bogdanov, Rybicki, & Grind- 
lay 2006, for examples of pulse profiles from non-magnetic hy- 
drogen atmosphere hot spots). We refer to the region at magnetic 
colatitude 0-10° as the primary magnetic polar cap and the op- 
posite pole as the secondary cap. The classification scheme (for 
isotropically-emitting hot caps) is defined in Beloborodov (2002): 
(I) only the primary cap is visible, and the pulse profile is purely 
sinusoidal with a single peak, (II) the secondary cap is seen around 
pulse minimum due to relativistic light-bending, which reduces the 
strength of the modulation, (111) the primary cap is not seen during 
a segment of the rotation, and (IV) both spots are seen at all phases 
and thus there is no modulation. Relativistic light-bending allows 
6 < 1 15° to be visible; thus, e.g., for (20°, 80°), one pole will not be 
seen during a portion of the rotation. 

Several important features are evident from a comparison of 
magnetic atmosphere emission to that of isotropic emission. The 
(St, 0t)-angular-dependence of the radiation (or beam pattern) man- 
ifests as a narrow "pencil-beam" along the direction of the magnetic 
field and a broad "fan-beam" at intermediate angles (see Pavlov 
et al. 1994; Lloyd 2003, for beam patterns and spectra at vari- 
ous ®B and (pt). As discussed in Pavlov et al. (1994), the pencil- 
beam is the result of the lower opacity at flj. < (is/fifl)''^, where 
Eg = heB/m^c = 11.6(i?/10'^ G) keV is the electron cyclotron 
energy; the width of the pencil-beam is thus ~ (ii/iis)''*, and the 




Figure 2. Light curves for different geometries (a, f): class I with (20°, 30°), 
class II with (30° ,60°), class III with (60°,80°), and class IV with (10°,80°). 
a is the angle between the spin and magnetic axes, and ( is the angle be- 
tween the spin axis and the line-of-sight. The four classes are defined in 
Beloborodov (2002). The solid lines are for the magnetic model described 
in the text [dashed line is for (50°, 50°)], while the dotted lines are analytic 
light curves (scaled arbitrarily in amplitude) for isotropic emission from two 
antipodal hot spots (see Beloborodov 2002). 



radiation is more strongly beamed at higher magnetic fields. For 
our case with B ~ 4 x 10'^ G, the width is ~ 6°. This narrow 
beam is seen in the (50°,50°)-light curve plotted in Figure|2l which 
is the only instance shown that has the observer's line-of-sight ex- 
actly crossing the magnetic cap and coinciding with the peak of 
the isotropic emission. Also evident is the fan-beam (most obvi- 
ous in the light curves of classes II and III), which occur on either 
side of the magnetic cap and can increase the number of light-curve 
peaks. In addition, the anisotropic beam pattern (combined with the 
surface temperature variation) can produce an apparent phase shift 
compared to isotropic emission and modulation when an isotropic 
beam pattern suggests none (c.f. class IV). 

Figure|3]shows the pulse fractions, PF = (Cmax-C'n,i„)/(Cn,ax + 
Cmin), where C is the count spectrum, as a function of energy for 
various geometries. The energy-dependence of the light curves and 
pulse fractions is due to the energy-dependent beam patterns and 
the surface temperature and magnetic field variations. As noted in 
Zavlin et al. (1995), the pulse fraction is lower at low energies since 
the beam pattern is more isotropic. Phase shifts between the peak 
of the low and high-energy light curves can also occur. 



4 COMPARISON TO RX J1856.5-3754 

TM07 discovered 7 s pulsations in the 2006 October 24 obser- 
vations of RX J1856.5-3754 made by the EPIC-pn and MOS 
cameras on XMM-Newton. These observations cover the energy 
range 0.15-1.2 keV. The pn light curve, with a pulse fraction of 
1.6%±0.2%, is shown in Figure|4]and found to have no significant 
difference when divided into soft and hard energies; the MOS light 
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Figure 3. Pulse fractions as a function of energy for different geometries 
(a,0 = (20° ,30°), (30° ,60°), (60°,80°), (10°,80°), (50° ,50°). a is tlie angle 
between the spin and magnetic axes, and f is the angle between the spin 
axis and the line-of-sight. 



Figure 4. Energy-integrated (0.15-1.2 keV) light curves of 
RX J1856. 5-3754. Histograms are the pn (top) and MOS (bottom) 
observations (see Fig. 1 of TM07). Solid and dotted lines are the models 
with (3°,40°) and (7°, 60°), respectively, along with the fit deviations [i.e., 
(data-model)/(T] in the corresponding lower panels. Two rotation periods 
are shown for claiity. 



curve is also shown in Figure|4]and found to be statistically consis- 
tent with no pulsations (TM07). 

As a result of this discovery, TM07 analyzed previous XMM- 
Newton observations and found light curves that are consistent be- 
tween the different observations; therefore these data were added 
together to obtain better statistics. Furthermore, when divided into 
soft (0.15-0.26 keV) and hard (0.26-1.2 keV) energy bands, the 
summed light curves demonstrate some energy-dependence, i.e., 
pulse fractions of 1.17% ± 0.08% (0.15-1.2 keV), 0.88% ± 0.11% 
(0.15-0.26 keV), and 1.5% ± 0.11% (0.26-1.2 keV; see TM07). 
The light curves are shown in Figure[5] 

Several characteristics of the observed light curves suggest 
possible values of a and ^: (1) the w 1% amplitude of the pulsa- 
tions, (2) a single peak (or two peaks close in phase, as may be 
the case for the pn light curve in Figure |4ll per rotation, and (3) 
no significant energy-dependence for the single observation (Fig- 
ure |4ll and small energy-dependence for the combined data (Fig- 
ure[5)- These characteristics imply the secondary cap is not a major 
contributor to the brightness variation (i.e., not class III) and small 
values of a or Large la - f | are also implied since otherwise 
the pulsations would be strong with a narrow peak when the cap 
becomes aligned with the direction to the observer (a similar argu- 
ment is made by Braje & Romani 2002 from the non-detection of 
RX J1856.5-3754 in the radio). 

We compare the light curves generated from our model first to 
the single (2006 October 24) XMM-Newton observation and then 
to the summed data. To avoid the computational tedious task of 
computing light curves for a continuous distribution in the angles a 
and (f, we select specific values of {a, i^), compute the light curve, 
and fit the model light curves to the observational data. The phase 
and normalization of the model light curves are arbitrary and are 
varied to minimize We restrict ourselves to integer values of 
(a, () since the uncertainties in our model (see H07 and Section lSTt 
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Figure 5. Light curves of RX J1856. 5-3754 for different energy ranges: 
0.15-1.2 keV (top), 0.15-0.26 keV (middle), and 0.26-1.2 keV (bottom). 
Histograms are the light curves obtained by summing together all the XMM- 
Newton observations (see Fig. 4 of TM07). The solid lines are the model 
with (2°, 30°), along with the fit deviations [i.e., (data-model)/cr] in the cor- 
responding lower panels, while the dashed lines are the model with (7°, 50°). 
Two rotation periods are shown for clarity. 
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Table 2. Fits to 2006 October 24 Light Curves 



Q,(a) 




^2/dof 


(deg) 


(deg) 


pn*) 




3 


20 


11.8/8 


8.12/9 


3 


30 


11.4/8 


8.36/9 


3 


40 


10.4/8 


6.72/9 


9 


50 


9.94/8 


9.92/9 


7 


60 


8.31/8 


9.25/9 


2 


70 


11.2/8 


7.48/9 



Notes: 

(a) a and f are interchangeable. 

(b) 10 phase bins - 2 fit parameters (phase and normalization) = 8 dof 

(c) The phase alignment from the pn fit is used. 

Table 3. Fits to Combined Light Curves 



(deg) (deg) 0.15-1.2 keV""' 0.15-0.26 keV^ 0.26-1.2 keV^ 



2 20 4.88/8 12.2/9 9.42/9 

2 25 5.42/8 9.95/9 7.07/9 

2 30 4.43/8 10.4/9 3.88/9 

2 35 5.10/8 9.36/9 4.27/9 

2 40 8.34/8 9.85/9 6.25/9 

4 45 6.93/8 16.7/9 9.49/9 



Notes: 

(a) a and f are interchangeable. 

(b) 10 phase bins — 2 fit parameters (phase and normalization) = 8 dof 

(c) The phase alignment from the 0.15-1.2 keV fit is used. 

and small variations between the data (c.f. Figures |4] and O do not 
warrant more detailed fits. 

We obtain good fits to the 2006 October 24 observation with a 
and ^ of ~ 2 - 9° and ~ 20 - 70°. Table[2]shows some of our fit re- 
sults, where the pn data is fit with a model light curve and then the 
MOS data is fit using the same model light curve and phase align- 
ment. Figure|4]illustrates two sample fits. These fits confine the al- 
lowed parameter space of (a, ^). Furthermore, the light curves from 
this observation do not show significant differences when divided 
into soft and hard energies, whereas the model light curves with 
^ > 50° show noticeably different hard-band pulsations. Therefore, 
the energy-(in)dependence of the data argues for ^ < 50°. 

Finally, we fit the (energy-dependent) light curves obtained 
by adding together the various XMM-Newton observations. Table[3] 
shows some of the results, where the total (0.15-L2 keV) data is 
fit with a model light curve and then the soft (0.15-0.26 keV) and 
hard-band (0.26- L2 keV) data is fit using the same phase align- 
ment. The more stringent constraints on (a, <f) are due to the better 
statistics and energy-dependence of the combined data. The model 
with (a, ()=(2°, 30°) is shown in Figure|5] 

As an illustration of our fit results, we show in Figure |6] the 
ranges allowed by the Icr uncertainty in the measured pulse frac- 
tions and the pulse fractions of the model light curves as a function 
of a and f . Only model light curves whose pulse fractions lie in 
the shaded regions can fit (within the uncertainties) the observed 
light curves. However, matching the pulse fraction is a necessary 
but not sufficient condition. As discussed above, the pulse shape 
and energy-dependence also serve as constraints. For example, it 
appears from Figure [6] that (7°, 50°) should fit the combined light 
curve. But the shape of the model light curve for this case is too 




a (deg) 

Figure 6. Pulse fractions as a function of f (or a) for different a (or ^). The 
horizontal lines indicate the pulse fractions from TM07 using data from the 
2006 October 24 observation and the combined observations; the shaded 
areas indicate the lo" uncertainty region. 



broad and shows large differences when divided into soft and hard- 
energy bands (see Figure [Sj. 



5 DISCUSSION 

In previous work (Ho et al. 2007), we successfully fit the multi- 
wavelength spectrum of RX J 1856.5-3754 with our model of ra- 
diation from a magnetic neutron star containing a partially-ionized 
hydrogen atmosphere. With the recent discovery of rotational mod- 
ulation of the X-ray emission (Tiengo & Mereghetti 2007), we use 
the light curves predicted by our model (which includes relatively 
simple surface distributions of the magnetic field and temperature) 
to constrain the geometry of RX J1856. 5-3754. In particular, we 
find angles of < 6° and » 20 - 45°: one is the angle a between the 
rotation and magnetic axes and the other is the angle ( between the 
rotation axis and the direction to the observer. Since the light curves 
with (a, if) are equivalent to those with (^, a), this implies that either 
the rotation and magnetic axes are closely aligned or we are essen- 
tially seeing down the spin axis of the neutron star (see also Braje & 
Romani 2002). We also note the following: (1) The model used here 
and described in Ho et al. (2007) assumes the high X-ray energy 
emission arises from a condensed iron surface below the hydrogen 
atmosphere. The emission properties of the condensed surface have 
been studied by Brinkmann (1980); Turolla, Zane, & Drake (2004); 
Perez-Azon'n, Miralles, & Pons (2005); van Adelsberg et al. (2005). 
(2) Though it would present an observational challenge, our model 
predicts optical pulsations with a singly-peaked sinusoid and pulse 
fractions of a; 0.1 - 0.5%. (3) Other members of the neutron star 
population RX J1856.5-3754 is thought to belong to show vary- 
ing degrees of pulsations, some with a single peak (per rotation) 
and others with double peaks. Analyses of these sources will be 
presented in future work. 
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5.1 Changes in Surface Distributions of B and T^g 

We point out an uncertainty in our model that is difficult to quan- 
tify but will affect the light curves and thus the quality of the fits 
and the angles {a, f ) inferred. We assume dipolar-like distributions 
of the magnetic field and temperature over the surface of the neu- 
tron star (see Table |T]l. If the size of each magnetic colatitudinal 
region is varied, the strength of the pulsations will change. For the 
geometries that best fit the observations, a +5° change in the size 
of the the region described by fl = 4 x 10'^ G, 0b = 60°, and 
Tcs = 5 X 10^ K can result in a < 0.002 change in the pulse frac- 
tion, which is on the order of the observational errors. In addition, 
if the distributions are described by an offset dipole or if quadrupo- 
lar components are added on top of the dipole, then the shape and 
strength of the pulsations will, in general, be different. 

5.2 RX J1856.5-3754 as a Normal Radio Pulsar? 

One set of our solutions (a < 6°) indicates a very close alignment 
between the rotation and magnetic axes. Observational evidence 
from studies of radio pulsars may also indicate the alignment of 
the rotation and magnetic axes on a timescale of ~ 10' y (Lyne & 
Manchester 1988; Tauris & Manchester 1998). For example, Tau- 
ris & Manchester (1998; see also Zhang, Jiang, & Mei 2003, and 
references therein) find the observed distribution of a peaks around 
30 - 50° ; after correcting for beaming, they find a decreases with 
time, with an average a x 40° for pulsars younger than 10^"' y. 
Thus both our solutions {a < 6° and a x 2Q - 45°) fit within this 
scenario, given the ~ 5 x 10' y age of RX J1856. 5-3754 (Kaplan, 
van Kerkwijk, & Anderson 2002; Walter & Lattimer 2002), 

In spite of its non-detection at radio frequencies (Brazier 
& Johnston 1999), suppose RX J1856.5-3754 is a normal ra- 
dio pulsar. The size of the radio beam p is estimated to be p a 
5 - 6°{PI\ s)-'/2 = 2° (see, e.g., Rankin 1993; Gould 1994). 
Our results allow the radio beam to miss our line-of-sight, with 
the (o- a 20 - 45°, f < 6°)-solution having a easier chance 
for this to occur. In addition, if RX J1856. 5-3754 is losing ro- 
tational energy by magnetic dipole radiation, the spindown rate 
is P = 10^'' s s-'(B/10'2 Gf(PI\ s)-' = 5 X 10^'' s s"', which 
is well below the P < 1.9 x 10"'- s s"' obtained by Tiengo & 
Mereghetti (2007). This spindown rate (and ~ 5x 10' y age) implies 
that RX JI856.5-3754 must have been born at roughly its cuiTent 
spin period if no significant change in P (or B) occurred over its 
lifetime. Also, an Ho' nebula was found around RX JI856.5-3754 
and has an estimated luminosity of fbow = 8 x 10^^- ergs s"' (van 
Kerkwijk & Kulkarni 2001b); the dipole spindown luminosity of 
RX J1856. 5-3754 is insufficient in powering this relativistic wind. 
However, it is important to note in the above considerations that, 
given its spin period and magnetic field, RX J1856. 5-3754 falls 
below the death line for radio pulsars; the death line being the re- 
gion in f - S-space (or P - P) where pulsars no longer emit in the 
radio (see, e.g., Sturrock 1971; Chen & Ruderman 1993; Harding, 
Muslimov, & Zhang 2002). 

5.3 X-ray Polarization 

In the presence of magnetic fields typical of isolated neutron stars 
such as RX JI856.5-3754 (B> 10'- G), radiation propagates in 
two polarization modes: the ordinary mode is polarized parallel 
to the k-B plane, while the extraordinary mode is polarized per- 
pendicular to the k-B plane (see, e.g., Meszaros 1992). Because 
the extraordinary mode opacity is greatly suppressed compared to 



the ordinary mode opacity [by a factor ~ (E/Eb)^', Lodenquai et 
al. 1974; Meszaros 1992], the extraordinary mode photons escape 
from deeper, hotter layers of the atmosphere than the ordinary mode 
photons, and the emergent radiation is linearly polarized to a high 
degree (as high as 100%; Gnedin & Sunyaev 1974; Meszaros et 
al. 1988; Pavlov & Zavlin 2000). Measurements of X-ray polariza- 
tion, particularly when phase-resolved and measured in different 
energy bands, could provide unique constraints on the magnetic 
field strength and geometry and the compactness of the neutron 
star (Meszaros et al. 1988; Pavlov & Zavlin 2000; Heyl, Shaviv, & 
Lloyd 2003; Lai & Ho 2003). As noted previously, light curves of 
the total flux with (a,() are equivalent to those with (^,cr). Thus 
the results presented in this work could be satisfied equally with 
an interchange of these two angles. However, the light curves of 
the polarized light for the two sets of angles are distinctly different 
(Pavlov & Zavlin 2000), with the ordinary mode spectra showing 
greater phase variability since the outer atmosphere layers are more 
sensitive to magnetic field orientation (Ventura et al. 1993). For our 
case, the polarization is essentially 100% in one polarization, and 
the position angle (in the plane of the sky) of the linear polarization 
would undergo small changes as the star rotates for the case with 
(2°, 30°), while the position angle would rotate 360° for the case 
with (30°, 2°). Propagation effects, such as those discussed in Heyl 
& Shaviv (2000, 2002); Heyl et al. (2003); Lai & Ho (2003), would 
need to be taken into account, but sources like RX JI856.5-3754 
could be studied by possible future X-ray polarization instruments, 
such as Constellation-X, XEUS, and the Extreme Physics Explorer 
(Bellazzini et al. 2006; Elvis 2006; Jahoda et al. 2007). 
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